Higher-Loop Structural Properties of the (5 Function in Asymptotically Free 

Vectorial Gauge Theories 
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We investigate some higher-loop structural properties of the /3 function in asymptotically free 
vectorial gauge theories. Our main focus is on theories with fermion contents that lead to an 
infrared (IR) zero in /3. We present analytic and numerical calculations of the value of the gauge 
coupling where ft reaches a minimum, the value of /3 at this minimum, and the slope of /3 at the IR 
zero, at two-, three-, and four-loop order. The slope of /3 at the IR zero is relevant for estimates 
of a dilaton mass in quasiconformal gauge theories. Some inequalities are derived concerning the 
dependence of the above quantities on loop order. A general inequality is derived concerning the 
dependence of the shift of the IR zero of /3, from the n-loop to the (n + l)-loop order, on the sign 
of the (n + l)-loop coefficient in /3. Some results are also given for gauge theories with M = 1 
supersymmetry. 



I. INTRODUCTION 



The evolution of an asymptotically free gauge theory 
from high Euclidean momentum scales u in the deep ul- 
traviolet (UV) to small scales in the infrared is of fun- 
damental field-theoretic interest. This evolution is de- 
scribed by the f3 function of the theory. Following the 
pioneering calculations of the j3 function at one- loop [1] 
and two-loop Q order, this function was subsequently 
calculated to three-loop @ and four-loop |4j order in the 
modified minimal [j| subtraction (MS) scheme @. The 
anomalous dimension of the (gauge-invariant) fermion bi- 
linear operator, j m , has also been calculated up to four- 
loop order in this scheme Q- 

Here we consider the UV to IR evolution of an asymp- 
totically free vectorial gauge theory with gauge group G 
and Nf massless fermions transforming according to a 
representation R of G 8]. An interesting property of 
this type of theory is that, for sufficiently large Nf, the 
two-loop /3 function has an IR zero @, If Nf is near 
to the maximum allowed by the property of asymptotic 
freedom, then this IR zero occurs at a small value, but, 
as Nf decreases, it increases to stronger coupling. This 
motivates the calculation of the IR zero of (3 at higher- 
loop order 10]. Calculations of this IR zero, and the 
associated anomalous dimension of the (gauge-invariant) 
fermion bilinear, j m , have recently been done to four- 
loop order for an asymptotically free vectorial gauge the- 
ory with gauge group G and Nf fermions in an arbi- 
trary representation R, with explicit results for R equal 
to the fundamental, adjoint, and symmetric and anti- 
symmetric rank-2 tensor representations flD . Il2j . A cor- 
responding analysis was carried out for an asymptotically 
free vectorial gauge theory with M = 1 supersymmetry 
in [13| . Although the terms in the j3 function at three- 



and higher-loop order, and the terms in 7 m at two- and 
higher-loop order are dependent on the scheme used for 
regularization and renormalization of the theory, these 
higher-loop calculations are valuable because they give a 
quantitative measure of the accuracy and stability of the 
lowest-order calculations of am and j m . A study of the 
effect of scheme transformations on results for am was 
performed in [l4j |. 

In this paper we will present calculations at the n-loop 
level, where n = 2, 3, 4, of several important quanti- 
ties that provide a detailed description of the UV to IR 
evolution of a theory with an IR zero in its /? function. 
Our general results apply for an arbitrary (non-Abelian) 
gauge group G. We denote the running gauge coupling at 
a scale « as gijj), and define a(fi) = g(^) 2 / (47r). (The fx 
argument will often be suppressed in the notation.) The 
loop order to which a quantity is calculated is indicated 
explicitly via the subscript n£, standing for n-loop, so 
that the n-loop ft function and its IR zero are denoted 
P n e and am,nt ■ Given the asymptotic freedom of the the- 
ory, the UV to IR evolution, as described by (3 n i, occurs 
in the interval 



< a(p) < am,, 



(1.1) 
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In addition to am. n i, the three structural properties of 
Pni that we study are (i) the value of a where j3 n t reaches 
its minimum in the interval (|1.1|) . denoted a m ,ni, (ii) the 
minimum value of j3 n i on this interval, (j3 n e)min, and 
(iii) the slope of f3 ne at a IR , n i, denoted dj3 n e/da\ aiR nf . 
The importance of the first two quantities for the UV to 
IR evolution of the theory is clear. One would like to 
know where the rate of running, (3 = da/dt, has max- 
imum magnitude, as a function of a, and hence, as a 
function of /i. Further, one is interested in what this 
maximum magnitude in the rate of running, i.e., (since 
/3 < 0), the minimum value of j3 is in the interval I a . 
The third quantity, the slope of the j3 function at am , is 
of interest because it describes how rapidly /3 approaches 
zero as a approaches am. A knowledge of this slope 
is also valuable because it is relevant for estimates of 
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a dilaton mass in gauge theories that exhibit approxi- 
mate scale-invariance associated with an IR zero of (3 at 
a value, am, that is sufficiently large that this approxi- 
mate dilatation symmetry is broken by the formation of 
a fcrmion condensate. For each of these structural quan- 
tities, one would like to see how higher-loop calculations 
compare with the two-loop computation. As part of our 
work, we derive some inequalities concerning the relative 
values of each of these quantities at the two- and three- 
loop order. 

We also generalize some results that were obtained in 
PH concerning am.ni- In [HI it was shown that in a the- 
ory with a given G, R, and Nf for which the two-loop (3 
function f3^i has an IR zero, the three-loop zero satisfies 
the inequality am 31 < am^i in the minimal subtrac- 



theory applicable in this region, one thus integrates out 
these now-massive fermions, and the (3 function changes 
to that of a pure gauge theory, which does not have any 
perturbative IR zero. Hence, in this case the infrared 
zero of (3 is an approximate, but not exact, fixed point of 
the renormalization group. 



tion (MS) scheme used there [15j. The reduction in the 
value of the IR zero going from two-loop to three-loop 
level is typically substantial; for example, for G = SU(2) 
and Nf — 8, ampt — 1-26, while am^t — 0.688, and 
for G = SU(3) and N f = 12, a m ,2t = 0.754, while 
a /-R,3f = 0.435. A natural question that arises from 
the analysis in [ill ] is how general this inequality is and, 
specifically, whether it also holds for other schemes. We 
address and answer this question here. We prove that 
for an asymptotically free theory with a given G, R, 
and Nf for which (3n has an IR zero, the inequality 
am y sz < am,2£ holds in any scheme which has the prop- 
erty that the sign of the three-loop coefficient in (3 is 
opposite to that of the one-loop coefficient for Nf G /, 
which thus preserves for fag the existence of an IR zero 
that was true of fiu>. This preservation of the two-loop 
IR zero in j3 is physically desirable, since /?2£ is scheme- 
independent, so if it exhibits an IR zero, then a reason- 
able scheme should maintain the existence of this zero at 
higher-loop level. More generally, we will derive a result 
that shows how am.ni shifts, upward or downward, to 
a iR,(n+i)e, when it is calculated to the next higher-loop 
order. 

For a given gauge group G, the infrared properties 
of the theory depend on the fermion representation R 
and the number of fermions, Nf. For a sufficiently large 
number, Nf, of fermions in a given representation (as 
bounded above by the requirement of asymptotic free- 
dom), the IR zero in (3 occurs at a relatively small value 
of a and the theory evolves from the UV to the IR with- 
out any spontaneous chiral symmetry breaking (S^SB). 
In this case, the IR zero of (3 is an exact infrared fixed 
point of the renormalization group. Thus, the infrared 
behavior of the theory exhibits scale-invariance (actu- 
ally conformal invariance [16j ) in a non- Abelian Coulomb 
phase. For small Nf, as the theory evolves from the 
UV to the IR, and the reference scale n decreases be- 
low a scale which may be denoted A, the gauge inter- 
action becomes strong enough to confine and produce 
bilinear fcrmion condensates, with the associated spon- 
taneous chiral symmetry breaking and dynamical gener- 
ation of fermion masses of order A. As /x decreases be- 
low A, and one constructs the effective low-energy field 



If Nf is only slightly less than the critical value Nf :Cr 
for spontaneous chiral symmetry breaking, so that am is 
only slightly greater than the critical value, a cr (depend- 
ing on G and R) for fermion condensation, then the UV 
to IR evolution exhibits approximate scale (dilatation) 
invariance for an extended logarithmic interval, because 
as a(/i) increases toward am, while less than a cr , j3 ap- 
proaches zero, i.e., the rate of change of a(/i) as a func- 
tion of n approaches zero. Thus, a(ji) is large, of 0(1), 
but slowly running ("walking"). This is quite different 
from the behavior of a s (/x) in quantum chromodynam- 
ics (QCD). This approximate scale invariance at strong 
coupling plays an important role in models with dynam- 
ical electroweak symmetry breaking [T3, EH , and occurs 
naturally in models with an approximate infrared fixed 
point [18|. Since am ~ O(l) and 7 m is a power se- 
ries in a, there is an enhancement of 7 m in such models, 
which, in turn, is useful for generating sufficiently large 
Standard Model fermion masses. Approximate calcula- 
tions of hadron masses and related quantities have been 
performed using continuum field theoretic methods for 
these theories . Recently, an intensive effort has been 
made using lattice methods to study the properties of 
SU(-/V C ) gauge theories with various fermion contents, in 
particular, theories that exhibit quasi-scale-invariant be- 
havior associated witn an exact or approximate IR zero 
of the respective j3 functions. For example, for SU(3) 
with fermions in the fundamental representation, mea- 
surements of 7 m have been reported in [2(j. In theories 
where N c , R, and Nf are such that am is only slightly 
greater than a cr , so this approximate scale invariance as- 
sociated with an IR zero of f3 at strong coupling holds, the 
spontaneous breaking of this symmetry by the formation 
of a fermion condensate, may lead to a light state which 
is an approximate Nambu-Goldstone boson (NGB), the 
dilaton [2l| (see also [l||). The mass of the dilaton de- 
pends on several quantities, including the effective value 
of the j3 function at the relevant scale fj, ~ A where the 
S^SB takes place. The desire to study the quasi-scale- 
invariant behavior of such a theory is an important moti- 
vation for obtaining more detailed information about the 
structure of the (3 function, as contained in the structural 
quantities (i)-(iii) discussed above. 
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II. BETA FUNCTION 
A. General 

The UV to IR evolution of the theory is described by 
the function 



Pa ~ dt ' 



(2.1) 



where t = In /z. This has the series expansion 

oo oo 

-2a^b,,a l , (2.2) 



= -2aJ2 b ea e 



where t denotes the number of loops involved in the cal- 
culation of be, a = 5 2 /(167r 2 ) = a/(47r), and 

(2.3) 



* (An) 1 

As noted above, the one- and two-loop coefficients b\ 
and &2j which are scheme-independent, were calculated 
in P, Q (see Appendix A) . The be for t > 3 are scheme- 
dependent; in the commonly used MS scheme, the be 
have been calculated up to four- loop order 0, E[. For 
analytical purposes it is more convenient to deal with 
the be , since they are free of factors of An. However, for 
numerical purposes, it is usually more convenient to use 
the be, since, as is evident from Table I of [ill ], the range 
of values of be is smaller than the range for the be- We 
will use both interchangeably. We denote the function 
calculated to n-loop order as 



t+i 



(2.4) 



c=i 



Some explicit examples of four-loop functions are given 
in Appendix B. 

With our sign conventions, the restriction to an asymp- 
totically free theory means that b\ > 0. This is equivalent 
to the condition 



Nf < N 



where \22L]2 



N 



f,bU 



f,blz 



IIC A 



AT 



f 



(2.5) 



(2.6) 



(biz stands for be zero). For the fundamental, adjoint, 
and symmetric and antisymmetric rank-2 tensor repre- 
sentations of G = SU(iV c ), the upper bound (12.5[) al- 
lows the following ranges of Nf. (i) Nf < (ll/2)iV c 
for fundamental, (ii) Nf < 11/ A for adjoint, (hi) Nf < 
IIN C /[2(N C ± 2)] for symmetric (antisymmetric) rank- 
2 tensor. In the case of a sufficiently large represen- 
tation R, this upper bound may forbid even the value 
Nf — 1. For example, for the rank-3 symmetric ten- 
sor representation of S\J(N C ), the upper bound is Nt < 



11 A^/ [(A c + 2)(N C + 3)], and the right-hand side of 
this bound is larger than 1 only for N c (analytically 
continued to nonnegative real numbers) in the interval 
3- V3 < N f <3 + V^, i.e., 1.268 < N f < 4.732 (to 
the indicated floating-point accuracy). Hence, if N c is 
equal to 2, 3, or 4, the bound allows only the single value 
Nf = 1, and if N c > 5, then the bound does not allow 
any nonzero (integer) value of Nf. For G = SU(2), with 
a representation labeled by the integer or half- integer j, 
the inequality (I2.5P is 



N f < 



33 



2j(j + l)(2j + l) 



for G = SU(2) . (2.7) 



This bound is: (i) N f < 11 if j = 1/2; (ii) N f < 11/ A if 
j = 1; (hi) N f < 11/10 if j = 3/2. The right-hand side 
of (12. 7[) decreases through 1 as j (continued to real num- 
bers) increases through 1.562, so that the upper bound 
(12. 7[) does not allow a nonzero number of fermions in a 
representation of SU(2) with j > 2 [H. 

To analyze the zeros of the n-loop function, n e, aside 
from the double zero at a = 0, one extracts the overall 
factor of —2a 2 and calculates the zeros of the reduced 
(r) polynomial 



0n£ 

2^ 



n 



a ' 



(2. 



or equivalently, be a . As is clear from Eq. 

the zeros of n e away from the origin depend only on 
7i — l ratios of coefficients, which can be taken as be/b n 
for I = l,...,7i— 1. Although Eq. (|2.8[) is an algebraic 
equation of degree n — 1, with n — I roots, only one of 
these is physically relevant as the IR zero of n e- We 
denote this as am tn e- In analyzing how the n-loop 
function describes the UV to IR evolution of the theory, 
we will focus on the interval (jl.lj) . 

To investigate how aiR,ni changes when one calculates 
it to higher-loop order, it is useful to characterize the full 
set of zeros of n i- In general, if one has a polynomial 
of degree to, P m (z) — ^2™ =Q k s z s , and one denotes the 
set of to roots of the equation P m (z) — as {zi, z m j , 
then the discriminant of this equation is defined as [251 ] 



(2.9) 



Since A m is a symmetric polynomial in the roots of 
the equation P m (z) = 0, the symmetric function theo- 
rem implies that it can be expressed as a polynomial in 
the coefficients of P m (z) [26]. We will sometimes indicate 
this dependence explicitly, writing A m (reo) K m ). The 
discriminant A m is a homogeneous polynomial of degree 
rn(m — 1) in the roots {z{\. For our present purpose, to 
analyze the zeros of n e away from the origin, given by 
the roots of Eq. (|2.8|) . of degree to = n — 1, we will thus 
use the discriminant A n _i(&i, b%, b n ), or equivalently, 
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A„_i(6i, 6 2 , —j b n ). Note that, because of the homogene- 
ity properties, 

An^ih, b 2 , K) = (47r)-(" +1 »"- 2 ) A„_i(6i, 62, 6„) . 

(2.10) 

Some further details on discriminants are given in Ap- 
pendix C. 

Although we focus on the behavior of (3 n £ in the phys- 
ical interval (jl.l[) . in characterizing the zeros of f3 n £, we 
will make use of some formal mathematical properties 
of f3 n i as an abstract function of a. For large \a\, since 
Pni ~ — 2&„a™ +1 , it follows that if a is large and posi- 
tive, then sgn(/3„f) = — sgn(6„), while for large negative 
a, sgn(/3„^) = sgn((-l)™6„). Thus, sgn(/3 n£ ) for large 
positive a is equal to (— l)" +1 sgn(/3„f ) for large negative 
a. Since f3 n £ is negative in the vicinity of the origin, it 
follows that for (both even and odd) n > 2, 

If b n < 0, then f3 n £ has at least one zero at a 
positive real value of a . (2-11) 

Furthermore, again because (3 ~ — 26„a" +1 for large \a\, 
a consequence is that for n > 2, 

If n is odd and 6„ < or n is even and 6 n > 0, then 



j3 n ii has at least one zero at a negative real 



value of a 



(2.12) 



Of course, the behavior of f3„i at negative values of a is 
not directly physical, and the behavior at large positive a 
is beyond the range of validity of the perturbative calcu- 
lation, but these mathematical properties will be useful 
in characterizing the total set of zeros of (3 n £ at higher- 
loop order. 

Given that Nf G I, so that P21 has an IR zero, we can 
track how this zero changes as the loop order n increases. 
One general result is as follows. As Nf /* Nf t },i z at the 
upper end of the interval /, a.iB. >n i — >• 0. This is a result of 
the fact that in this limit, b\ — >• 0, so that pni,r reduces 
to ai^™ =2 c/ _2 , which has a root at a = 0. Starting 
at the (n — 2)-loop level and tracking the physical IR 
zero at three- and higher-loop order, one can infer that 
generically a n ,e is the root of Pni,r that moves toward 
zero in this limit Nf /* Nfj,iz- 

Because /3 n £ is a polynomial in a(p) and hence a con- 
tinuous function, and because /3 n £ = at the two ends 
of the interval (jl.ip . at a = and a = am tn i, and is 
negative for small (positive) a, it follows that j3 reaches 
a minimum in this interval (ll.l[) . This occurs at a point 
where df3 n e/da = 0, which we label u m .ni (where the 
subscript m stands for "minimum /3 in / a "), and we de- 
note 



a—a m n i 

From Eq. (|2.4p . one calculates d(3 n e/da 



(2.13) 



(47r) 1 d/3 n i/da^ with the result 

d/3 n i 



da 



2^+ l)b e a £ = -2^2(1 + l)b t <^ . (2.14) 



The equation for the critical points, where df3 n i/da = 0, 
is thus an algebraic equation of degree n, with n formal 
roots, one of which is a = 0. Assuming that 62 < 0, 
i.e., Nf G /, so that the two-loop /3 function has an 
IR zero (and also, in higher-loop calculations, that the 
scheme preserves the existence of this IR zero), it follows 
that, among the remaining n — 1 roots, one is real and 
positive and yields the minimum value of f3 n £ for a in 
the relevant interval (jl.ip . and this root is the above- 
mentioned OL mn £. 

Given that f3 has an IR zero at czir and is analytic at 
this point, one may expand it in a Taylor series about 
am. This involves the slope of the /3 function at am. 
For compact notation, we denote 



8> =^ 
PlR - da 



and, for the n-loop quantities, 

at _ df3 n e 



J IR,nl 



da 



(2.15) 



(2.16) 



With [3 (am) = 0, the expansion of (3(a) for a near to 
a IR is 

/3 = (3' m (a - am) + o((a - a IR ) 2 ) . (2.17) 

Here we have written this expansion for the full (3 func- 
tion; a corresponding equation applies for /3 n £. 



B. IR Zero of /3 at the Two-Loop Level 

We next review some background on the two-loop j3 
function that is relevant for present work. The two-loop 
f3 function is P21 = —2a 2 (61 + b%a). This has an IR zero 
at 



am,: 



bi _ Anbi 
b 2 b 2 



(2.18) 



which is physical if and only if 6 2 < 0. The coefficient 62 is 
a linear, monotonically decreasing function of Nf, which 
is positive for zero and small Nf and passes through zero, 
reversing sign, as Nf increases through Nf^2z, where 



N 



fMz 



2T f (5C A + 3Cf) 



(2.19) 



For arbitrary G and i?, Nf t bu > Nf^z, as is proved by 
the fact that 



_ 3C A (7C A + llC f ) 
" N f , b 2z - 4T/(5Ca + 3C/) > 



(2.20) 
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Hence, there is always an interval I of Nf values for which 
the two-loop (3 function has an IR zero, namely 

/ : AW < N f < N fMz . (2.21) 

For example, in the case of fermions in the fundamental 
representation, denoted □, 

34W 3 1 1 N 

6U\ C lliV. i£R = n (2 22) 

13iV c 2 - 3 1 2 v ' 

so that, for N c = 2, the interval I is 5.55 < Nf < 11; 
for N c = 3, I is 8.05 < Nf < 16.5; and as N c -> oo, / 
approaches 34JV C /13 < Nf < HJV c /2. 

Since we are primarily interested in studying the IR 
zero of P and since the presence or absence of an IR zero 
of the two-loop {3 function, P21, is a scheme- independent 
property, we focus on Nf G /, where this IR zero of /?2£ is 
present. A general result is that for a given gauge group 
G and fermion representation R and Nf S /, am,it is 
a monotonically decreasing function of Nf. As Nf de- 
creases from Nf^iz, aiR y 2Z increases from 0. As Nf 
decreases through a value labeled Nf lCr , am increases 
through a critical value, a cr ~ 0(1), where fermion con- 
densation takes place. Thus, 

Nf = N fiCr ^ am = a cr . (2.23) 

The value of Nf jCr is of fundamental importance in the 
study of a non-Abelian gauge theory, since it separates 
two different regimes of IR behavior, viz., an IR confor- 
mal phase with no SxSB for Nf >cr < Nf and an IR phase 
with SxSB for Nf < Nf )Cr . As Nf approaches Nf^iz a * 
the lower end of the interval /, am,2i becomes too large 
for Eq. (|2~T8l) to be reliable. 

Because of the strong-coupling nature of the physics 
at an approximate IR fixed point with am ~ 0(1) j 
there are significant higher-order corrections to results 
obtained from the two-loop B function, which motivated 
the calculation of the location of the IR zero in /?, and 
the resultant value of 7 m evaluated at this IR zero, to 
higher-loop order for a general G, R, and Nf [TTL [T3| . 

C. Function and Dilaton Mass 

Here we focus on a theory in which the IR zero of (3, 
am, is slightly greater than a cr , so that, in the UV to IR 
flow, there is an extended interval in t = In \i over which 
a(fi) is approaching am from below, but is still less than 
a cr . In this interval, a(p) ~ 0(1), but /3 is small, and 
hence the theory is approximately scale-invariant. As /1 
decreases through A, a(/i) increases through a(A) = a cri 
the fermion condensate forms, and the fermions gain dy- 
namical masses, this approximate scale invariance is bro- 
ken spontaneously. In terms of the (symmetric) energy- 
momentum tensor 6^, the dilatation current is = 
d» u x v , and one has d^D" = [fi / {4a)}G a ^G afw , where 
G^„ is the field-strength tensor for the gauge field. When 



taking matrix elements, the deviation of this divergence 
d^D^ from zero, i.e., the nonconservation of the dilata- 
tion current, thus arises from two sources, namely the 
facts that /? is not exactly equal to zero and the nonzero 
value of the matrix element of G^„G aAI!y , defined appro- 
priately at the scale A. An analysis of the matrix element 
of T>^ between the vacuum and the dilaton state |x(p)), 
in conjunction with a dimensional estimate of the gluon 
matrix element, and the Taylor series expansion (|2.17p 
evaluated with /1 ~ A yields the resulting estimate for 
the dilaton mass m x [2l| 

m 2 x ~ p' IR (a m - a cr )A 2 . (2.24) 

In terms of n-loop level quantities, the right-hand side of 
Eq. (|2.24p /3' IR ni i a iRM — a cr )A 2 . The importance of 
the slope at am, (3' IR , and the n-loop calculation of this 
slope, f3' IR nl , in estimating a dilaton mass in a quasicon- 
formal theory is evident from Eq. (I2.24p . As is the case 
with am,,ni, because of the strong-coupling nature of the 
physics, it is valuable to compute higher-loop corrections 
to the two-loop result, f3' IR 2i ■ Below, we will present 
two- and higher-loop analytic and numerical calculations 
of (3' IR nt . Other effects on m x have been discussed in the 
literature [2lj . including the effect of dynamical fermion 
mass generation associated with the spontaneous chiral 
symmetry breaking as /1 descends through the value A. 
Owing to this and other nonperturbative effects on m x , 
we restrict ourselves here to presenting one input to this 
calculation, namely f3' IR nl , for which we can give definite 
analytic and numerical results. 



D. IR Zero of (3 at the Three-Loop Level 

Let us assume that Nf € /, so that P21 has an IR zero. 
Here we analyze how this IR changes as one calculates 
the B function to three- loop order, extending our results 
in [1 11 ] to (an infinite set of) schemes more general than 
the MS scheme used in that paper. Since the existence 
of of the IR zero in the two-loop j3 function is a scheme- 
independent property of the theory, it is reasonable to 
restrict to schemes that preserve this IR zero of j3 at the 
three-loop level. We first determine a condition for this 
to hold. 

The three-loop B function is /33£ = — 2a 2 /?3£, r , so, aside 
from the double zero at a — (the UV fixed point), 
vanishes at the two roots of the factor Pzi,r = b\+ 62a + 
63 a = 0, namely, 



at/3z,3e,± = 




= ^(-6 2 ±VA 2 (6i, b 2 ,6 3 ) ) -(2-25) 

where A 2 (6i, 62,63) = b 2 — 46^3 . The analysis of the 
IR zero of requires an consideration of the sign of 
A2 (61, 62, 63). The condition that j3u have an IR zero 
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requires, in particular, that its two zeros away from the 
origin be real, i.e., that A 2 {b\,b 2 , b 3 ) > 0. For a given G, 
R, and Nf £ /, so that b\ and b 2 are fixed, this condition 
amounts to an upper bound on b 3 , namely b 3 < 6|/(4fei). 
Now, b 2 — > at the lower end of the interval /, so that, 
insofar as one considers the analytic continuation of Nf 
from positive integers to positive real numbers, the above 
bound generically requires that 63 < for Nf £ I. This is 
also required if R = □, and one studies the theory in the 
limit N c — > 00 and Nf — > 00 with r = Nf/N c fixed, since 
in this case there are discrete pairs of values (N c ,Nf) 
that enable one to approach arbitrarily close to the lower 
end of the interval I at r = 34/13 where b 2 —> 0. In 
order to preserve the existence of the two-loop IR zero at 
the three-loop level, one is thus motivated to restrict to 
schemes in which 63 < for Nf £ /, and we will do so 
here. (The marginal case 63 = is not generic, since 63 
varies as a function of iV c and Nf, so we will not consider 
it further.) 

Before proceeding, it is worthwhile to recall how the 
property 63 < for Nf £ / arises in the MS scheme. 
In this scheme, 63 is a quadratic function of Nf with 
positive coefficients of the N 2 term and the term in- 
dependent of Nf. This coefficient 63 vanishes, with 
sign reversal, at two values of Nf, denoted iV/,f,3 Z ,- and 
Nf,b3z,+i given as Eq. (3.16) in [ll|, with 63 < for 
Nf, b 3z - <N f < N mZi+ (and 63 > for Nf < N mz - 
and Nf > N ftb3Zl+ . In [ll| it was shown that in this 
scheme, for all of the representations considered there, 
namely, the fundamental (□), adjoint, and rank-2 sym- 
metric ( l I l ) and antisymmetric (H) tensor representa- 
tions, Nf t b3 Z - < Nf : b2z and Nf lb 3z,2 > Nf t bu, so that 
b 3 < for all Nt £ I. For example, for fermions in 
the R = □ representation, (i) for N c — 2, Nf t b3z,i = 
3.99 < N f . b2z = 5.55 and N f . b3z , 2 = 27.6 > N fMz = 11; 
(ii) for N c = 3, N ftb3Ztl = 5.84 < N ftb2z = 8.05, 
and iV/,b3 2 ,2 = 40.6 > N fMz = 16.5; (iii) as N c -> 
00, Nf' b 3 Z ,i -> 1.911iV c while N ftb2z -> 2.6157V C and 
Nf, b3z ,2 -> 13.348iV c , while N m z ^ 5.5N C . In Table 
|TT]we list the values of A2(6i, 62, ^3) with b 3 calculated in 
the MS scheme, for the illustrative cases N c — 2, 3, 4 
and Nf in the respective / intervals. Since b 3 < for 
iV/ € / in this MS scheme, it follows that all of the 
entries in this table have A 2 > 0. 

Given that 63 < for Nf € I, Eq. (|2.25p can be rewrit- 
ten as a = (2iT/\b 3 \)(-\b 2 \Ty / bl + 46i|6 3 |). The solution 
with a — sign in front of the square root is negative and 
hence unphysical; the other is positive and is am t3 £ , i.e., 



&IR,3 



2tt 
N 



46i|6 3 



(2.26) 



In [Uj it was shown that in the MS scheme, for all 
Nf £ I, am, 3 e < aiR i2 £. Here we demonstrate that 
this result holds more generally than just in the MS 
scheme. We prove that for arbitrary gauge group G, 
fermion representation R, and Nf £ I, in any scheme 
in which 63 < for Nf £ I (which is thus guaranteed 



to preserve the IR zero present at the two- loop level), it 
follows that ajji 3 ( < ajR 2 £. To prove this, we consider 
the difference 



2tt 
I&2&3 



26i|6 3 



±h\b 3 



(2.27) 

The expression in square brackets is positive if and only 
if 



(2&1I&3 



bUbl + ^M) >o 



(2.28) 



This difference is equal to the nonnegative quantity 
(2bib 3 ) 2 , which proves the inequality. Note that, since 
b\ is nonzero for asymptotic freedom, this difference van- 
ishes if and only if 63 = 0, in which case am i3 i = am t2 £. 
We have therefore proved that 

ai R ^£ < a m . 2 £ if 63 < for Nf £ I. (2.29) 

As noted above, am^ 2 £ is a monotonically decreasing 
function of Nf £ I. With 63 < for Nf £ I, this mono- 
tonicity property is also true of am t3 £. As Nf increases 
from Nf^ b2z to Nf^ b \ z in the interval /, aiR t u decreases 
from 



ciiRM = ^\lih &tNf = Nf, b 2z (2.30) 

to zero as Nf /*■ Nf lb iz at the upper end of this interval, 
vanishing like 



OilRM 



Anbi 



1 



+ 0(6?) 



(2.31) 



as Nf /• N fMz and 61 -> 0. 



E. IR Zero of /3 at the Four-Loop Level 



The four- loop j3 function is fiu = — 2a 2 (3^^, so ftu 
has three zeros away from the origin, at the roots of the 
cubic equation 



^ 4 £,r = h + b 2 a + b 3 a 2 + b A a A = , 



(2.32) 



(where fi n £. r was given in Eq. (|2.8|l ). These zeros were 
analyzed for the MS scheme in [IJ, [l3] . Here we extend 
this analysis to a more general class of schemes that have 
b 3 < for Nf £ I, and hence maintain at the three-loop 
level the IR zero of the scheme-independent two-loop (3 
function. 

The nature of the roots of Eq. ()2.32|) is determined by 
the sign of the discriminant A3 (61, b 2 , b 3 , 64), or equiva- 
lently, 

A 3 = A 3 (b 1 ,b 2 ,b 3 ,b 4 ) = b 2 2 b 2 3 - 27b 2 1 bl-4(b 1 bl + b 4 bl) 



I8M2M4 



(2.33) 
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The following properties of A3 are relevant here: (i) if 
A 3 > 0, then all of the roots of Eq. (|2.32l) are real; (ii) if 
A3 < 0, then Eq. (|2.32[) has one real root and a complex- 
conjugate pair of roots; (hi) if A3 = 0, then at least two 
of the roots of Eq. (|2 . 32[) coincide. Given the scheme- 
independent properties 61 > and b 2 < (i.e., Nf € I), 
and provided that 63 < for Nf £ I, we can write this 
discriminant as 

A 3 (6i, b 2 ,b 3 , 64) = b 2 b 2 -27b 2 b 2 +4(b 1 \b 3 \ 3 + b 4 \b 2 \ 3 ) 
+ 186 1 |6 2 ||6 3 |6 4 ■ (2.34) 

If 64 were zero, then the zeros of f3u would coincide with 
those of /?3£, and the property that these are all real is in 
accord with the reduction 

A 3 (h ,b 2 ,b 3 ,0) = b 2 3 (b 2 2 + 4h 1 63 1 ) = b 2 3 A 2 (h , b 2 , 63) , 

(2.35) 

which is positive. 

Now consider nonzero 64. First, assume that the 
scheme has the property that 64 > 0. Then we can write 
Eq. ([2321 as 

h - \b 2 \a- \b 3 \a 2 + ha 3 = if 64 > . (2.36) 

From an application of the Descartes theorem on roots 
of algebraic equations, it follows that there are at most 
two (real) positive roots of this equation and at most one 
negative root. Moreover, from Eq. (12.12)) we can deduce 
that in this case with 64 > 0, in addition to the double 
zero at a = 0, fin has a zero at a negative value of a, 
so the upper bound on negative zeros from the Descartes 
theorem is saturated. Furthermore, since f3u is negative 
at large positive a, there are then two possibilities: cither 
the two remaining zeros of Eq. (|2.32l) arc a complex- 
conjugate pair, or else they are both real and positive. 

If, on the other hand, the scheme is such that 64 < 0, 
then we can write Eq. (|2.32j) as 

h - \h\ot- \b 3 \a 2 - \h\a 3 = if 64 < . (2.37) 

From a similar application of the Descartes theorem, we 
infer that there is at most one positive real root and at 
most two negative real roots of Eq. (12.371) . From Eq. 
(|2.11[) we deduce that fta has a zero at a positive real 
value of a. Depending on I&4J, the other two roots of Eq. 
(|2.37p may be real and negative or may form a complex- 
conjugate pair. 

Combining the information from both the Descartes 
theorem and the discriminant A3 , we derive the following 
conclusions about the roots of Eq. (|2.32p and hence the 
zeros of j3u aside from the double zero at a = 0. As 
before, we assume that b 2 < (i.e., Nf 6 7) so that fce 
has an IR zero, and also that the scheme is such that 
b 3 < for Nf € I, guaranteeing that this IR zero is 
maintained at the three-loop level. Then, 

1. If 64 > and A 3 > 0, then Eq. (|2~32"[) has one 

negative and two positive real roots, 



2. If 64 > and A 3 < 0, then Eq. (|2.32l) has one neg- 
ative root and a complex-conjugate pair of roots, 

3. If 64 < and A 3 > 0, then Eq. $Z3Z§ has one 
positive root and two negative roots, 

4. If 64 < and A 3 < 0, then Eq. ([232]) has one 
positive root and a complex-conjugate pair of roots. 

For a particular pair (JV C , Nf), the marginal case A 3 = 
might occur, and would mean that two of the roots of 
Eq. (|2.32[) are degenerate. Since this equation is a cubic, 
it would follow that all of the roots are real. If A3 = 
and 64 > 0, then Eq. (|2 . 32[) has one negative root and 
a positive root with multiplicity 2, while if 64 < 0, then 
(|2.32|) has one positive root and a negative root with 
multiplicity 2. 

It is reasonable to avoid schemes that lead to the out- 
come (2) above, with no real positive root of Eq. (|2.32[) . 
since these fail to preserve the IR zero of the scheme- 
independent two-loop (3 function. Although the positiv- 
ity of A3 is not a necessary condition for this preserving 
of the IR zero, it is a sufficient condition. We thus investi- 
gate the conditions under which A3 is positive. As shown 
via Eq. I|2.35p . if 64 = 0, then A 3 > 0. By continuity for 
small 1 64 1, A3 remains positive, and there is only a small 
shift in the two zeros that were present in fi 3 g, together 
with the appearance of a new zero. Since the highest- 
degree term in A 3 (bi,b 2 ,b 3 ,b4) involving |&4|, namely 
—21b\b\ is negative-definite, it follows that, other things 
being equal, for sufficiently large |£>4 ) , A 3 (b\,b 2 ,b 3l b^) 
will decrease through zero and become negative. The 
two 64 values at which A3(&i, b 2 , 63, 64) = are 

\b 2 \(2b 2 + % 1 \b 3 \)±2(b 2 + 3b 1 \b 3 \) 3 / 2 

(MA3*,± = ^ 

(2.38) 

Therefore, a sufficient condition for a scheme to be such 
that /?4£ preserves the IR zero that is present in (3 2 £ and 
Pu for Nf e I is 

(h) A 3z- < 64 < (Ma3*,+ . (2.39) 

Note that at the lower end of the interval /, where b 2 — s- 0, 
the interval ([2.391) reduces to the upper bound |&4| < 
2|6 3 | 3 / 2 /(276 1 ) 1 /2. 

For reference, in the MS scheme, 64 is a cubic polyno- 
mial in Nf and is positive for Nf <E I for N c — 2,3 but 
is negative in part of I for higher values of N c (see Table 
I of [ll|, where N_ c is denoted N), In Table [U we list the 
values of A 3 (bi,b 2 ,b 3 ,bi) with 63 and 64 calculated in 
the MS scheme, for the illustrative values N c = 2, 3, 4 
and values of Nf in the respective I intervals. For all 
of the four-loop entries in Table U A 3 (bi,b 2 ,b 3 ,l>4) > 0, 
as is evident from the values listed explicitly in Table [Til 
so these entries correspond to the case (1) in the list of 
possibilities for 64 and A 3 given above. 

Rather than calculating aiR,ni directly from @ n i, a 
different approach is to use (3 n i to compute Pade approx- 
imants and then calculate zeros of these approximants. 
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As before, since one is interested in the zeros away from 
the origin, one extracts the factor —2a 2 in Eq. (|2.4p and 
analyzes the polynomial fi n t,r m Eq. (|2.8|) . of degree n— 1 
in a, depending on the n coefficients be, £ — 1, .., n. From 
this, one can construct a set of \p, q] Pade approximants, 
i.e., rational functions, each with a numerator polyno- 
mial of degree p and a denominator polynomial of degree 
q in a, of the form (J2^ =0 Pj^ j )/(J2k=o1kZ k ). Without 
loss of generality, one can divide numerator and denomi- 
nator by qg, so that, after redefinition of the coefficients, 
the [p, q] Pade approximant to Eq. (I2.8P is 



i + E 



(2.40) 



depending on the p + q + 1 coefficients pj with j = 0, ...,p 
and qfc with fc = 1, .., q. These coefficients arc determined 
by matching the Taylor series expansion of \p, q] in a with 
the n coefficients be, £ = 1, n, so that p + q = n — 1. 
Thus, from the four- loop beta function factor fiu,r, one 
can construct two relevant approximants with p + q = 3, 
namely the [1,2 ] and [2,1] Pade approximants [27j]. We 
did this in [ll| and calculated the resultant unique IR 
zero from the [1,2] approximant and the relevant IR zero 
from [2,1], denoted otin : At, [1,2} and a iR,Ae, [2,1}, respec- 
tively. These were found to be close to the directly cal- 
culated IR zero, otiR,&t.. From the known results for (3 n e 
with n = 2, 3, 4, one can make estimates of f3§e by 
various methods, but since these only contain exact in- 
formation up to the n = 4 loop level, we will not pursue 
this direction here. 

As the n = 4 special case of the result discussed above, 
o-ir,u decreases to zero as Nf f~ Nf^u- For the R = □ 
and for a given Nf g / where it is reliably calculable, 
c-ir,ai is slightly larger than aj R ^e, but the difference, 
o-iR,At — ctiR,3£ is sufficiently small that aiR^e is smaller 
than aiR^i- For higher fermion representations and Nf 
values where the IR zero is reliably calculable (i.e., not 
too close to the lower end of the interval /), the difference 
ctiRfli — otiR^i is again smaller in magnitude then the 
difference olir^i — &iR,3e but may have either sign. Thus, 
where aiR^e is reliably calculable, it is smaller than 0,2,1- 
The finding that the fractional change in the location of 
the IR zero of (3 is reduced at higher-loop order agrees 
with the general expectation that calculating a quantity 
to higher order in perturbation theory should give a more 
stable and accurate result. 

The scheme-dependence of aiR tTl e for n > 3 can be 
studied by carrying out scheme transformations, recal- 
culating a' IR ni in the new scheme, and comparing with 
o-iR,ni- This study was carried out in [14| . To be accept- 
able, a scheme transformation must satisfy a number of 
necessary conditions, such as mapping a a positive real a 
to a positive real a' and vice versa. Although these condi- 
tions can be satisfied easily in the vicinity of the ultravio- 
let fixed point of an asymptotically free theory at a — 0, 
they are nontrivial and constitute significant restrictions 
on scheme transformations at an infrared fixed point [l4j . 
For example, the scheme transformation a — tanha', 



with inverse a' = (l/2)ln[(l +a)/(l — a)], is acceptable 
for small a, in the vicinity of the UV fixed point of an 
asymptotically free gauge theory, but is not acceptable in 
the vicinity of an IR fixed point at a = aiR ~ O(l), since 
a can approach 1 from below, in which case a' diverges, 
and a can exceed 1, in which case a' is complex. 

Scheme-dependence of higher-loop calculations is 
present not just in calculations of an IR zero of f3 n e 
at three- and higher-loop level, but also in higher-loop 
perturbative QCD calculations. The fact that the MS 
scheme is a reasonable one has been demonstrated, e.g., 
by the excellent fit that has been obtained to experi- 
mental data for a s (p) with fi 2 = Q 2 using this scheme 
[28| . There has been much work on optimized schemes 
for higher-order QCD calculations [2^. However, we note 
that two of the simplest scheme transformations that one 
might apply for QCD are not generally acceptable at an 
IR zero of j3 with ajR ~ O(l). These are the scheme 
transformations denoted S2 and S3 in (T3 |. which are 
constructed to render the leading scheme-dependent co- 
efficient in the new scheme, b' 3 , equal to zero. They are 
acceptable at the UV zero of /3 and hence in perturbative 
QCD applications, but are not, in general, acceptable in 
the vicinity of an IR zero with aiR ~ 0(1) because they 
can map a real positive a in the MS scheme to a negative 
or complex coupling in the transformed scheme, as was 
shown in 1141. 



F. Shift of IR Zero at (n + l)-loop Level 

Here we derive a result on the direction of the shift in 
the IR zero of the /3 function when one increases the order 
of calculation of /3 from the n-loop level to the (n+l)-loop 
level, where n > 2. We assume, as before, that the theory 
is asymptotically free and that 62 < (i.e., Nf g I), so 
that there is an IR zero of /3 at the two-loop level. We 
assume that the scheme-dependent coefficients be with 
£ = 3, n + 1 are such that they preserve the existence 
of the IR zero of (3 at higher-loop level [30]. We focus 
here on values of a close to aiR, n e, where d/3 n e/da > 0. 
Expanding /3 n e in a Taylor series expansion around a = 
otiR,nl, with the abbreviation P' IR nl = d(3 ne /da\ aiR ni 
defined above, we write the general Eq. (|2.17|) explicitly 
in terms of n-loop quantities as 

fine = &iR, n i ( a - a m,ni) + 0((a - a /fi ,„£) 2 ) . (2.41) 

Now let us calculate f3 to the next-higher- loop order, i.e., 
fi(n+i)ii and solve for the zero, ct IR ^ n+1 ^, which cor- 
responds to aiR tTl £ (among the n — 1 zeros of P( n +i)i 
away from the origin). To determine whether am ( n +i)t 
is larg er or smaller them cxjr i.e., whether there is a 
shift to the right or left, consider the difference 



P{n+i)i ~ fine — —2b n+ ia 



n+2 



(2.42) 



In a scheme in which b n+ \ > 0, this difference, eval- 
uated at a = ajR.ni, is negative, so, given that 
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d(3 n z/da\a IR nl > 0, to compensate for this, the zero 
shifts to the right, whereas if b n+ i < 0, the difference 
is positive, so the zero shifts to the left. That is, 

If b n+ i > , then a IR ^ n+1 y > ai R . n t 

If b n+l < , then a IR ^ n+1)e < a IR . n g . (2.43) 

(In a scheme with 6 n +i = 0, obviously aiRt n +i)e, — 
ctiR,n£)- The application of this general result (|2.43p is 
evident in the specific calculations in at the three- 
loop and four-loop levels. 

III. fi FUNCTION STRUCTURE 

At high scales in the UV, the fi function is dominated 
by the leading quadratic term, fi ~ — 2b\a 2 + 0(a 3 ). 
The calculation of the IR zero of fi n i is important for 
investigating the UV to IR evolution of the theory. But, 
as discussed in the introduction, for a more detailed study 
of this evolution, one needs not just the value of the IR 
zero, cti Rt ni, but the full curve of fi n g for ol G I a . Here we 
present calculations of three quantities that give further 
information about this curve, including (i) the value of 
a where fi n i reaches its minimum for a G I a , a mt ni, (ii) 
the minimum value of fi n e for a G I a , (fi n e)min', and 
the slope fi' IR nl at the IR zero of fi, as defined in Eq. 
(I2.16[) . The relevance of the third quantity to estimates 
of the dilaton mass in a quasiconformal gauge theory has 
been noted above. Our calculations are performed at the 
n = 2, n = 3, and (n = 4)-loop level. 



A. Two-Loop Level 

1. Position of Minimum in fiii for a 6 I a 

At the two-loop level, given that 62 < so that 
fill function has an IR zero, the derivative dfiiijda = 
—2a(2bi + 3&2o0 vanishes at a = a mt 2t, where 



cx>m.ii 



26i 
'362" 



87r6i 87rfei 



Sb 2 



3I&2 



Explicitly, 



8tt(11C a - 4T f N f ) 
3[4(504 + 3Cf)T f N f - 34C1] 



(3.1) 



(3.2) 



327r(llCU — ATfNf) 3 
81[MC 2 A ~4(5C A + 3C f )T f N f ] 



Note that 



,2£ — ^ a IR,2t 



3. Slope of f^2t at aiR^t 



(3-3) 



(3.4) 



The derivative dfoi/da evaluated at a = ai R ^.i is 



IR.11 



26| 



26? 



26? 

N 



2(11C A -4T / JV / ) S 



3[A(5C A + 3C f )T f N f -3AC 2 A ] 



(3.5) 



which is positive for Nf G /. 

As descriptors of the shape and structure of the (3 func- 
tion, the quantities a mM , (/3 ni ) min , and P' IR nl are inter- 
related. Thus, if one makes a rough, linear (lin.) approx- 
imation to the (3 function in the interval from a — a m , n e 
to a = aj R ^ n £, then this slope would be 



A« 



Ot-IR,; 



(3.6) 



For example, in the (n = 2)-loop case, substituting the 
values of (fiii) m in, olir.u, and a m .2ti this approximation 
yields 



A/3: 



2ldin. 



%2 



9\h 



(3.7) 



which exhibits the same dependence on the input coef- 
ficients bi and 62, with a somewhat smaller coefficient, 
8/9 rather than the coefficient 2 in the exact two-loop 
expression, (3' IR2 ei m Eq. (|3.5I) . 

In Tables HI IIII1 and IIVI we list numerical values of 
a m ,niy {Pni)min, and P' IR n( for fermions in the R = □ 
representation of SU(A^), for some illustrative cases of 
N c and, for each N c , values of Nf in the respective in- 
tervals /. As illustrations, we show plots of j3 n i in Fig. 
[T] for iV c = 2 and N f = 8 and in Fig. 1 for iV c = 3 
and Nf = 12 as functions of a. The results in the tables 
and figures are given for the quantities evaluated at the 
n = 2, n = 3, and n = 4 loop levels. The n = 3 and 
n = 4 loop results will be discussed further below. 



2. Minimum Value of P21 for ot G I a 

At a = a m ^£, @2t reaches its minimum physical value 
for a G I a , namely 



(Mr 



86? 32tt6? 



276i 



2761 



B. Three-Loop Level 

1. Position of Minimum in fizi for a £ I a 

Here we assume a scheme in which 63 ^ 0, since if one is 
working with a scheme in which 63 = 0, then (3^ = fin, so 
the analysis of the three-loop fi function reduces to that 
of the two-loop fi function discussed above. Furthermore, 
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FIG. 1: Plot of the n-loop function n i as a function of a for 
n = 2, 3, 4 and N c = 2, Nf = 8 with fermions in the fundamental 
representation. At a given value of a, the curves, from bottom to 
top, are for 02£, 04,1, and 0^£, respectively. See text for further 
details. 



for the reasons explained above, we restrict to schemes 
in which 63 < for Nf £ /. The derivative dfizi/da = 
— 2a(26i + 36 2 a + 463a 2 ) is zero at a = and at the two 
other points, 



a = 



1 

8h 



36 2 ± \ %i - 326163 



(3.8) 



This can be written as a 



7 r(2|6 3 |)- 1 (-3|6 2 | T 

y / 96f+ 3267|6 3 | ) . The critical point corresponding to the 
— sign in front of the square root is negative and hence 
unphysical, while the critical point corresponding to the 
+ sign in front of the square root is a m ^£, i- e -, 



&m,3£ 



3I&2 



32b 1 \b 3 



2|&3. 

A general inequality is 

We prove this by examining the difference 

7r r 

a m ,2t- a m,u = g ,, r 1 166i|6 3 |+96 
6|6 2 6 3 | L 



(3.9) 



(3.10) 



- 3|6 2 |^/962 + 326 1 |6 3 | j (3.11) 

The condition that this is positive is equivalent to the 
condition that the square of the polynomial term in the 



FIG. 2: Plot of the n-loop function n £ as a function of a for n = 
2, 3, 4 and the illustrative case N c = 3, Nj = 12 with fermions in 
the fundamental representation. At a given value of o, the curves, 
from bottom to top, are for 02i, 04i, and 0^, respectively. See 
text for further details. 



2. Minimum Value of 0zi for a 6 I a 

At a = ctj3 mini 3£, foe reaches its minimum value for 
a <E I a , namely 



(Mr 



64|6 3 



(l446i62|6 3 | + 1286 2 62 + 2764 



|6 2 |(96 2 + 326 1 |6 3 



3/2 



(3.12) 



Note that one can write (Pu)min in terms of the bi co- 
efficients by replacing each be in Eq. (|3. 16[) by the cor- 
responding bi and dividing the overall expression by 4tt. 
Since (fi n i)min < 0, it is convenient to deal with the 
magnitudes \(j3 n i)min\- We find the following general in- 
equality: for a given G, R, and Nf £ I, in a scheme that 
has 63 < and hence maintains the existence of the IR 
zero in /3 2 £, 



I (Aw) I < \{P2f)min\ ■ 
To prove this, we consider the difference 



(3.13) 
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I (/Mr 



\(Mr. 



1728^2 1 6 3 1 3 



20486f |6 3 | 3 + 27bj(U4b 1 b 2 2 \b 3 \ + 128bfb 



2^2 



27$ 



27|6 2 | 3 (96l 



326!|6 3 



\3/2 



(3,141 



The positivity of this difference is equivalent to the pos- 
itivity of the square of the polynomial terms in the nu- 
merator minus the square of the term in the numerator 
involving the radical. This difference of squares is equal 
to 

81926f |6 3 | 3 (5126f |6 3 | 3 +34026^&i|6 3 |+17286l6f 6§+7296l) . 

(3-15) 

This expression is manifestly positive-definite, which 
proves the inequality (|3.13|) . 

3. Slope of p 3t at a IR - ie 

The derivative of f3 3 i at a = am^i is 
1 



ft 



1R.M 



\b 3 ? 



4|&2|(&2 + &lM) 



(6 2 + 26 1 |6 3 |) 



4&1I63 



(3.16) 



That this is positive follows from the fact that the square 
of the term involving the square root minus the square 
of — (4foi I &2 1 1 4*3 1 + l^l 3 ) in the brackets is the manifestly 
positive quantity 46f &§(fo| + 46i|6 3 |). Owing to the ho- 
mogeneity properties, to express fi' IR 3e in terms of the bt 
coefficients, one simply replaces each bg in Eq. (|3.16[) by 
the corresponding bg. 

A general inequality is that for a given G, R, and Nf £ 
/, in a scheme with 63 < 0, which is thus guaranteed to 
maintain the existence of the IR zero in P21 at the three- 
loop level, 



Pirm < ft 



IR,2t 



(3.17) 



To prove this, we examine the difference 



ft 



IRSlt 



ft 



IR,U 



2b\b, 



2 U 2 



•4|62|(^ + 6i|6a|)-N(^ 



26i|6 3 |) 



46i|6 3 



(3.18) 



r 



The positivity of this difference is equivalent to the pos- 
itivity of the square of the polynomial terms in the nu- 
merator minus the square of the term in the numerator 
involving the square root, which is 



Abjbj + \2b\b\b\ 



156? 



24^61 1 63 



(3.19) 



This is manifestly positive-definite, which proves the in- 
equality (|3~T7l> . 

The shifts in the values of the IR zero, cxzR,nit the 
position of the minimum in (3 n i, the value of f3 n e at the 
minimum, and the slope of (3 n i at a = a n ^ are evident 
from Tables U QUI and El and Figs. ED and EJ 



C. Four-Loop Level 



in Tables HI IIII1 and II VI These four-loop structural results 
are also evident in Figs. [T]and[2j 

In addition to the results that we have proved above, 
we note some others here. As stated above, numerical re- 
sults for three- and four-loop structural quantities were 
calculated in the MS scheme. First, although the ra- 
tios a m ^i/ otiR,u and OL m ^ij am,M are not constants as 
functions of Nf, they do not differ very much from the 
two-loop ratio, which is a constant, namely, 2/3, as given 
in Eq. (13.41) . With fermions in the □ representation, for a 
given N c and Nf G /, a m ^ is slightly larger than a m ^, 
but still substantially smaller than a m ^.i^ just as is true 
of the corresonding am tn e quantities. Moreover, for a 
given G and loop order n, a m ^ is a monotonically de- 
creasing function of Nf G / and vanishes as Nf /• Nf^\ z 
and b\ — > 0. 



_The derivative d/3 3e /da = -2u(2b x + 3b 2 a + Ab 3 a 2 + 
564a 3 ) is zero at a — and at the three other points 
giyen by the zeros of the cubic equation 2bi + 3l>2a + 
Ab 3 a 2 + 564a 3 = 0. We have calculated these critical 
points, evaluated fta at its minimum physical value, and 
also evaluated the derivative dp^/da at a — am^g. We 
give the numerical results for a m ^£, (ft^min, and j3' IR u 



IV. SOME PROPERTIES OF -y m 

The anomalous dimension j m for the fermion bilinear 
t/jip describes the scaling properties of this operator and 
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can be expressed as a series in a or equivalently, a: 

oo oo 



1=1 



where q — ci/(Att) 1 is the £-loop series coefficient. The 
coefficient c\ is scheme-independent, while q for i > 2 
are scheme-dependent. The q coefficients have been cal- 
culated up to four-loop order in the MS scheme [zj- We 
list ci for £ — 1, 2, 3 in Appendix A. We denote the 
n-loop expression for "f m as a series in a, evaluated at 
the n-loop IR zero of j3, a — ot.iR >n i-> as Iirm- 

In [ll[ , we calculated ^iR.ni up to (n = 4)-loop order in 
the MS scheme. An important result was that we found 
a substantial reduction in jjr going from the two-loop 
to three-loop level for all of the fermion representations 
that were considered. The difference going from three- to 
four-loop level, "fiR t 3i — Jir,m> was found to be smaller 
and could be of either sign, depending on the represen- 
tation and value of Nf. The resultant ^ir^u was thus 
substantially smaller than jiR^e- 

One may investigate the reduction Jir^u < ^irm 
found in [1J| further. To do this for a given gauge 
group G and fermion representation R, we assume that 
Nf E I, so that the theory has an IR zero of /J-y, and, 
further, that the scheme is such that 63 < for Nf 6 /, 
so that this IR zero is guaranteed to be maintained at 
the three-loop level. We will use the resultant prop- 
erty that aiR^e < o*-ir,ii- Let us consider the difference 
Iirm - JiR,3t- This is given by 

lIR,2i — lIR.U = Ci(aiR,2£ — 0-IR,m) 



c 2 (a 



IR,U 



(4.2) 



The (scheme-independent) coefficient c\ is positive, so 
that, since ajR,2£ — oiir,m > 0, it follows that the first 
term on the right-hand side of Eq. (|4.2p is positive. The 



factor (aj R 2 e" a iR. 3/) m the second term is also positive. 
The coefficient C2 is scheme-dependent, so the analysis of 
this term necessarily involves a choice of scheme, as does 
the analysis of the third term. We next prove that in 
the MS scheme, C2 > for all of the representations 
considered in [ll|, so that this second term is positive. 

I 



To show this, we begin with the □ representation, for 
which 

_ (jV c 2 -l)(203A c 2 -9-20A e iV / ) 

\ C 2)fund,MS — 192A 2 ' 

The first factor in the numerator, N^ — 1, is obviously 
positive for all physical iV c . The second factor is positive 
for Nf < Nf tC 2 Z , where 



N 



203A, 2 



f,c2z 



20N C 



(4.4) 



This is larger than the upper bound on Nf from asymp- 
totic freedom, Nf^i z , as is clear from the difference 



N 



f,c2z 



N 



f,blz 



so that Cr, 



3(31A C 2 - 3) 
20N. 



> 



(4.5) 



l, 2 fund jjg > for Nf £ I (actually for all physi- 
cal Nf). This provides an anal ytic understanding of the 
numerical results in Table V of [ll| , which indicated that 
c 2 > for all N c and Nf considered there. 

We next consider the case of fermions in the adjoint 
representation, for which 



A c 2 (53 - 10Nf) 
24 



(4.6) 



This is positive for Nf < 53/10, which is larger than 
the upper bound on Nf for this representation from the 
requirement of asymptotic freedom, namely Nf^iz — 
11/4 = 2.75, so again, c 2 > for all Nf for the adjoint 
representation in the MS scheme. 

Finally, we consider the case of fermions in the sym- 
metric or antisymmetric rank-2 tensor representation, de- 
noted 5*2 and A2, respectively, with Young tableaux CD 
and Owing to the fact that the A2 representation 
of SU(2) is the singlet, it is understood that N c > 3 in 
this case. Since various formulas are similar for these 
two representations, with appropriate reversals of signs 
of certain terms, it is convenient to give them in a unified 
manner, with T2 referring to 5*2 and A2 together. We 
have 



(N c ± 2){N C T 1) 109A 2 ± 9A C - 18 - 10A C (A^ C ± 2)N f ] 



-2,T2,MS 



48 A 2 



(4.7) 



r 



where the upper (lower) sign applies for the S2 (A2) rep- 
resentation, respectively. In the numerator of this ex- 
pression, the factor ( A c ± 2) (N c =p 1) is obviously positive 
for the relevant values of N c , so one next examines the 



factor [109A 2 ± 9Ac - 18 - 10A c (Ac ± 2)N f }. This is 
positive for Nf < Nf tC 2T2z, where 



A 



f,c2T2z 



106A 2 ± 9A C - 18 
10A C (A C ± 2) 



(4.8) 
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As for the other representations, Nf :C 2T2z is larger than 
the respective upper bound on Nf from the requirement 
of asymptotic freedom, Nf t biz,T2, 



N f 



ll(N c 



«f M *,T2- 2{Nc±2) ■ 

This is proved by considering the difference 



(4.9) 



3(17iV c 2 ±3iV c -6) 

Nf }C 2T2z - Nf, b lz,T2 = 77777 tAT — ^ (4.1U) 



10JV c (iV c ± 2) 



This difference is positive for 



N r > 



34 



(4.11) 



i.e., 0.5124 for S2 and 0.6888 for A2, and hence for all 
physical N c . Therefore, this proves that c 2 > for all 
relevant Nf < Nt^iz and, in particular, for all Nf in the 
respective intervals / for these theories with fermions in 
the symmetric or antisymmetric rank-2 representation. 

We have thus proved that for the MS scheme, for all 
of the representations considered in (Tl| . the first two 
terms in the difference "fiR t2 £ — 1ir,u are both positive. 
We have also investigated the contribution of the third 
term. By analytic methods similar to those exhibited 
above, we find that this third term also makes a positive 
contribution to the difference in Eq. (|4.2I) . i.e., C3 < for 
Nf € I, in most, although not all, cases. For example, 
for G = SU(iVc) and fermions in the □ representation, 
c 3 < for all iV c up to N c — 15 and integer Nf values in 
the respective intervals /. This includes all of the cases 
of N r and Nf € / for which numerical results were given 
in [1 11 ] and thus gives an analytic understanding of those 
results. For N c = 16, the interval / is 42 < Nf < 87, 
and C3 < for all of these values of Nf except the lowest 
one, Nf = 42, where C3 > 0. Similar comments apply for 
larger N c . 



V. SUPERSYMMETRIC GAUGE THEORY 
A. IR Zeros of ft 

It is of interest to give some corresponding results on 
properties of the ft function and associated UV to IR 
evolution in an asymptotically free, M = 1 supersym- 
metric gauge theory with vectorial chiral supcrficld con- 
tent (f>i, <ti, i = l,...,Nf in the R, R representations, 
respectively. A number of exact results have been de- 
rived describing the infrared properties of the theory in 
[31], [32j]. Thus, one can compare findings from perturba- 
tive calculations with these exact results, and this was 
done in [l3j]. The ft function of the theory has the form 



da 
lit 



1=1 



a 



(5.1) 



where we use the subscript s, standing for supersymmet- 
ric, to avoid confusion with the corresponding quantities 
in the nonsupersymmetric theory, and bi, s = 6^ s /(47r)^. 
The beta function calculated to n-loop order is denoted 
ftne,s- For values of Nf for which ft n ^ s has an IR zero, 
we denote this as ain, n e,s- In addition to the scheme- 
independent coefficients b\^ s and &2,s! calculated in (33| 
and [3J], respectively, the three- loop coefficient, b^.i, has 
been calculated in [35| in the dimensional reduction (DR) 
scheme (36|. Calculations of azn, n ^ 3 and corresponding 
values of the anomalous dimension for the bilinear chi- 
ral superfield operator $<!> were given in [l3| up to the 
maximal order to which bg and the coefficients of the 
anomalous dimension had been calculated, namely the 
three-loop level. 

We recall that, since b\ = 3CA~2TfNf [331 ]. the upper 
bound on Nf for the theory to be asymptotically free is 



Nf < Nf lb u 



where 



N 



3CU 



/,61z,s 



2T 



f 



(5.2) 



(5.3) 



The two-loop ft function coefficient is 62, s = 6C\— 4(Ca+ 
2Cf)TfNf [34|, which decreases through positive values 
and passes through zero, reversing sign, as Nf increases 
through 



N 



f,b2z,s 



2T f (C A + 2Cf) 



(5.4) 



Since Nf^ 2 z,s < -^/.fciz.s, there is always an interval of 
values of Nf, namely 



N ftb2ZtS <N f < N f , bU 



(5.5) 



in which the two-loop ft function for this theory has an 
IR zero for Nf. The value of this two- loop IR zero is 



4tt6 m 

2tt(3C a - 2T f Nf) 
2{C A + 2Cf)T f Nf - 3C2 



(5.6) 



In particular, for chiral superfields in the □ and □ repre- 
sentations, N fMZjS = 3N C and N ftb2z , s = 37V c 3 /[2jV c 2 - 1], 
so 



37V, 3 



2N? - 1 



< N f < 



3N r 



for R = □ 



(5.7) 



In this case, the exact value of Nf at the lower end of 
the IR-conformal, non-Abelian Coulomb phase was de- 
termined in Ref. 1321 to be 



N 



/,cr,s 



3N C 
2 



for R = □ 



(5.8) 
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Here, since Nf t b2z,s > Nf tCr , s , the coefficient b 2 , s passes 
through zero and reverses sign in the interior of the non- 
Abelian Coulomb phase. Consequently, as was noted 
in [la ], for this case of chiral superfields in the □ and 
□ representations, one cannot study the IR zero of (3 2 e 
throughout the entirety of this phase. The generalization 
of iVj iCr s to higher representations has been given as 37] 



2T f 



(5.9) 



Assuming that j m saturates its upper bound of 1 in this 
supersymmetric theory as Nf \ Nf :Cr ^ s , Eq. (|5.8[) agrees 
with the result obtained via a closed-form solution for /3 
[3lj | (see also [11]). Note that the fact that the scheme 
used in [3l[ is different from the DR scheme does not 
affect this, since Nf >cr>s is a physical quantity. 

At the three- loop level, (3u,s = — ^oc 2 f3u,r, Sl where 
P n e,r is given by Eq. (|2.8I) with bi replaced by bg tS . One 
makes use of the result [35| 

b 3 , s = 21C 3 A + 4(-5C 2 A -13C A Cf+4Cj)T f Nf 



+ 4(C A + 6C f )TfN 



(5.10) 
H. One 



in the DR scheme. The three-loop IR zero of f} t 
was calculated and compared with am^e.s m 
can prove various inequalities similar to those that we 
have proved above for a non-supersymmetric gauge the- 
ory. We illustrate one of these, concerning the relative 
size of otiR^t,s and a/i? ; 3£ iS for chiral superfields in the □ 
and □ representations. We begin by noting that 63,3 in 



the DR scheme is a quadratic function of Nf which is pos- 
itive for small Nf, and, as Nf increases, passes through 
zero, becoming negative, at a value denoted Nf t b3z,-, s , 
reaches a minimum, and then passes through zero again 
at Nf — Nf_b3z,+,s, and is positive for larger Nf. In 
general, 



N 



f,b3z,±,s 



hC\ + l3C A C f 



fit. 



2T f (C A +6Cf) 



(5.11) 



where 



R s = AC\ + 4C%C f + \29C 2 A C) - 10AC A C 3 f + 16C f 



For the □ or □ representation, this reduces to 



4 

(5.12) 



21N? - 9N 2 - 2 ± JR S fund 



2N C (AN 2 - 3) 



where 



Rsjund = 105A^ C 8 - 126A^ C 6 - 3A^ C 4 



36iV; 



(5.14) 



(Note that R s j U nd is positive-definite, vanishing at eight 
complex values of Nf.) To prove that am^e.s < otiR,M,s 
for this case, it suffices to show that &3 )S < for Nf G I s , 
since then one can apply the same proof that we used for 
Eq. ([2~29|) . To show that b 3:S < for Nf G I s , we will 
demonstrate that Nf^3 Z ,-, s < Nf ib2ZiS and Nf t b3z,+,s > 
Nf,biz,s- First, for this case of R equal to the □ repre- 
sentation, we consider the difference 



J 



Nffi2z,s — Nf } b3z 



-187V C 6 + 21iV c 4 - 5AT? - 1 + (2N 2 - l)y/R s 



fund 



2N C (2N 2 - 1)(47V C 2 - 3) 



(5.15) 



r 



Although the polynomial term in the numerator of (|5.15[) 
is negative, it is smaller than the term involving the 
square root. To show this, we observe that the square 
of the term involving the square root minus the square 
of the polynomial term in the numerator is 24Nc 4 (N 2 + 
1)(N 2 - 1) 2 (47V C 2 - 3). This is positive for all physical N c , 
proving that Nf^2z,s > Nf^3z,~,s for this case. Next, we 
consider the difference 



N 



f,b3z,+,s 



N 



-3N? 



f,bU 



9N 2 -2 + y/R, 



s,fund 



2N C (4N 2 - 3) 

(5.16) 

Although the polynomial term in the numerator is nega- 
tive for physical N c , it is smaller than the square root, as 
is shown by the fact that the difference of the square of 
the square root term minus the square of the polynomial 
term is 247V 2 (TV 4 - 1)(4A C 2 - 3), which is positive. So 
we have proved that for this case, Nf t b3 Z ,-,s < Nf t b2z,s 



and Nf : b3z,+.s > ^V/,6iz,s- In turn, this proves that for 
this case with Nf chiral superfields in the □ and □ repre- 
sentations, in the DR scheme, 63^ < for Nf G I s , and 
hence 



OiIR,3Ls < &IR,2L. 



(5.17) 



This inequality follows by the same type of proof as the 
one that we gave for Eq. (|2 .29[) . 



B. Structural Properties of f3 3 

Because one has exact, nonpcrturbative results avail- 
able for this theory, we will be brief in our discussion of 
structural properties of the (3 function. The value of a 
where fi 2 £,s has zero slope and a minimum in the interval 
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FIG. 3: Plot of the n-loop function n i a for an SU(A r c ) gauge 
theory with M = 1 supersymmetry, as a function of a, for n = 2 
and n = 3 loops and N c = 2, Nf = 5 with chiral supcrficlds 
in the fundamental representation. The lower and upper curves 
correspond to @2£,s and fe^si respectively. See text for further 
details. 



FIG. 4: Plot of the n-loop (5 function f3 n £ g for an SU(iV c ) gauge 
theory with M = 1 supersymmetry, as a function of a, for n = 2 
and n = 3 loops and N c = 3, Nt = 7 with chiral superfields 
in the fundamental representation. The lower and upper curves 
correspond to /32i yS and /?3f lS , respectively. See text for further 
details. 



([1~T|) is given by Eq. d3TTJ) as 

8tt(3CU - 2TfNf) 



3[4(CU + 2C,)T f N f - 6C\] ' 
At this a, /?2£,s reaches its minimum physical value 



(5.18) 



32tt(3Ca ~ 2T f N f 



27[4(C A + 2C f )T f N f -6C A ] 



212 . (5.19) 



As in the nonsupersymmetric theory, a mt 2t. s — 
(2/3)ai R , 2 e,s- 

The derivative dj32i,s/da evaluated at a — aiR^e.s is 
given by the analogue of Eq. (I3.5[) . namely 



IR,2i,s 



2{iC A - 2T / A / ) 2 
4(C^ + 2C»T / JV> - 6C1 ' 



(5.20) 



which is positive for iV* € /. 

At the three-loop level, a m , 3 £, s , (/33<?)m»n, and P' IR3itS 
are given by Eqs. (|3T9]l . (|3TT2"j) . and ([3~T6l> with the re- 
placements bi — ¥ bi^ s . 

In Figs. [3] and |4] we show plots of the two- loop and 
three-loop (5 functions for this supersymmetric gauge the- 
ory with chiral superfields in the fundamental representa- 
tion and with the illustrative values N c = 2, Nt = 5 and 
N c = 3, Nf = 7, respectively. The three-loop j3 functions 
are calculated in the DR scheme. 



VI. CONCLUSIONS 

In this paper we have studied some higher-loop struc- 
tural properties of the j3 function in an asymptotically 
free vectorial gauge theory, focusing on the case where 
the theory has an IR zero in the f3 function. These struc- 
tural properties include the value of a where /3 reaches 
a minimum (i.e., a maximal magnitude, since /? < 
for a € I a), the value of /3 at this minimum, and the 
derivative d/3/da at the IR zero, calculated to the n- 
loop order. We have given results up to four loops in a 
non-supersymmetric gauge theory and up to three loops 
in a gauge theory with J\f = 1 supersymmetry. In an 
asymptotically free theory with an exact or approximate 
infrared zero in the j3 function, these structural quanti- 
ties provide further information about the running of a 
as a function of the reference scale, \x. The derivative 
of j3 at am is also of interest because it enters into es- 
timates of the dilaton mass in a quasiconformal gauge 
theory. A general inequality was proved concerning how 
the shift in the IR zero of /3 as one goes from the n-loop 
to the (n + l)-loop order depends on the sign of b n+ i. For 
schemes which have 63 < for Nf G I and which thus 
are guaranteed to preserve the existence of the IR zero 
in the (scheme-independent) P21 at the three-loop level, 
we have proved that o,ir,m < ceiR,2£, a mM < a m , 2 <?, 
■nun I < \(M min 

I, and 0' IRU < P' IR2 f Our re- 
sults further elucidate the ultraviolet to infrared evolu- 
tion of an asymptotically free vectorial gauge theory with 
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fermions. 



In the M S scheme 3 
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2857 
~5T 



C 3 A + +T f N f 



2C) 



205 



C A C f 



~tT Ca 



(TfN f f 



44 ,158 

Y c f + —C A 



(7.3) 



We have also used the four- loop coefficient, 64, calculated 
in the MS scheme in for our calculations. This coef- 
ficient 64 is a cubic polynomial in Nf. 



VII. APPENDIX A 



A. P Function Coefficients 



B. Coefficients for y„ 



For a vectorial gauge theory with gauge group G and 
Nf fermions in the representation R, the coefficients 61 
and 62 in the j3 function are [l[ 



We list here the c e for £ = 1, 2, 3: 



d = 6C f 



(7.4) 



h = -(llC A -4T f Nf) 



and @, [H 



(7.1) 



b 2 = -[34Cl-A(5C A + 3Cf)T f Nf] . (7.2) 



c 2 = 2C f 



3 97 10 

2 C/ + y CA -y T/iV/ 



(7.5) 



c 3 = 2C/ 



129 



129 



11413 
108 



Ci + Cf(TfN f ) ( - 46 + 48C(3)) - C A (T f N f ) + 48C(3)) - ^ (TyJV/) 



140 



(7.6) 

where £(s) = X^^Li n s i s the Riemann zeta function, and £(3) — 1.2020569.... We have also used the four-loop 
coefficient C4, calculated in the MS scheme in [7], for our calculations. 



VIII. APPENDIX B 

Here we give some illustrative explicit (3^ functions for 
various values of N c and Nf. The three- loop and four- 
loop coefficients are calculated in the MS scheme. These 
are written in the form 



and are listed both analytically and numerically (to the 
indicated floating-point accuracy). 



p u = -2h a 2 



1=2 1 



(8.1) 



N c = 2, N f 



a 

7T 



603 /a\ 2 
64 V7J 



-0.3183a 2 ( 1 - 0.7958a - 0.9546a 2 



9216 
0.3562a 3 



■13685!)- l!)sr,2N L ,(3)^aj3 



•2) 



N c = 2, N f = 9 



B 



169 /a 
~ ~32~W 



154445 / a\ 2 



9216 



(;)' 



-22506041 + 50531904C(3)\ /a\ 3 



1327104 



)(;)' 
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-0.2122a 2 (l - 1.6811a- 1.6980a 2 + 0.9292a : 



■3) 



iV c - 3, N f = 10 : /3 4 



13a 2 

67T 



37 /a 
1 ~ 26 W 



41351 /a\ 2 /-13418011+ 13331592C(3)\ /a\3 



3744 



(;)+( 



404352 



(=)' 



0.6897a 2 1 - 0.4530a - 1.1191a 2 + 0.2080a £ 



(8.4) 



A^ c = 3, Nf = 12 : p u 



3a_ 2 r 
' 2tt 



25 /a\ 6361 /a\ 2 / -140881 + 219192C(3)\ /a 
~6~ W ~ 288 W + V 3456 /W 



-0.4775a 2 (l - 1.3263a - 2.2379a 2 + 1.1441a £ 



(8.5) 



IX. APPENDIX C 



For m = 3, 



Consider the polynomial of degree to in z, P m (z) = 
S^Lo Ks2;S - As discussed in Section HH information on 
the nature of the roots of the equation P m (z) = is 
given by the discriminant A m defined in Eq. (|2.9[) . Since 
A m is a symmetric function of the roots (being propor- 
tional to the square of the Vandermonde polynomial of 
these roots) , the theorem on symmetric functions [26[ im- 
plies that A m can be expressed as a polynomial in the 
coefficients k s , s = 0, ...,m. We indicate this in the no- 
tation A m = A m (nQ, K m ). The discriminant A m is 
most conveniently calculated in terms of the Sylvester 
matrix of P(z) and dP(z)/dz, equivalent to the resultant 
matrix, denoted 5p,p', of dimension (2to— 1) x (2m — 1): 



S 



P 3 ,Pi 



so that 



A 3 (rto, «1,K2, «s) 



^3 


3k 3 





«'2 

2k 2 
3k 3 




K 2 
Kl 

2k 2 
3k 2 



«0 



2k 2 



(kik 2 ) 2 - 27(k K3' 

18kqKiK 2 K3 . 






Kl ) 



(9.3) 



2 - 4(k k 2 + k 3 kD 
(9.4) 



— /I \m(m-l)/2 -1 



A m = (-1 



det(5p,p/) . 



(9-1) 



Since we will use A m for m 
explicit expressions here: 



2 and to = 3, we list the 



A 2 (k , ki, k 2 ) = k 2 - 4k k 2 



(9.2) 
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TABLE I: Value of ct mn i a * the n-loop level with n = 2, 3, 4 
for an SU(A r c ) gauge theory with Nf fermions in the fundamental 
representation, with Nf £ /. As discussed in the text, a m „£ is the 
value at which the n-loop f3 function takes on its minimum value in 
the interval < a < ajR n i. Results are given for the illustrative 
values N c = 2, 3, 4. For comparison, we also list the IR zeros of /3 
calculated at n-loop level, ctjR n (, for n = 2, 3, 4, from Ref. [ill . 
For this and other tables, quantities evaluated at the n = 3 and 
n = 4 loop level are calculated in the MS scheme. 
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Otm,2l 


OL m ;it 
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10 
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1.47 


0.534 
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0.402 
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0.470 
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0.317 


0.337 
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14 
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TABLE II: Values of the discriminants A2(fci, 62, £3) and 
^3(61 , 62, 63, ^4) (see Eqs. 1 12. 96 and 1 12. lOt ) for the three-loop and 
four-loop IR zero equations, with B3 and B4 calculated in the MS 
scheme. Results are given for the illustrative values N c = 2, 3, 4. 
Notation ae-n means a X 10 _n . 
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TABLE III: Minimum value of the n-loop /3 function, /3 n £, denoted 
(flnl)min, in the interval < a < ajR^ n i relevant for the UV to 
IR evolution, calculated to n = 2, 3, 4 loop order for an SU(7V C ) 
theory with Nf fermions in the fundamental representation, with 
Nf 6 /. Values are given for N c = 2, 3, 4. Notation ac-n means 
a X 10-". 
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TABLE IV: Value of df3 n i/da. at n = 2, 3, 4 loop order for an 
SXJ(N C ) theory with Nf fermions in the fundamental representa- 
tion, with Nf £ /, evaluated at the IR zero calculated to this order, 
a lR,nl- We denote this here as j3' IR n( . 
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